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CMB Theory from Nucleosynthesis to Recombination 



1 Introduction 

These lecture notes comprise an introduction to the well-established physics and phenomenology of 
the cosmic microwave background (CMB) between big bang nucleosynthesis and recombination. We 
take our study through recombination since most of the temperature and polarization anisotropy 
observed in the CMB formed during the recombination epoch when free electrons became bound into 
hydrogen and helium. 

The other reason for considering only this restricted range from nucleosynthesis to recombination 
is that these notes are meant to complement the other lectures of the XIX Canary Island Winter 
School of Astrophysics. While they are self-contained and complete in and of themselves, they omit 
important topics that will be covered elsewhere in this volume: namely inflation (Sabino Matarrese), 
observations (Bruce Partridge), statistical analysis (Licia Verde), secondary anisotropy (Matthias 
Bartelmann), and non-Gaussianity (Enrique Martinez-Gonzalez). 

Furthermore the approach taken here of introducing only as much detail as necessary to build 
physical intuition is more suitable as a general overview rather than a rigorous treatment for training 
specialists. As such these notes complement the more formal lecture notes for the Trieste school 



which may anachronistically be viewed as a continuation of these notes in the same notation (Hu 



2003| ). 

The outline of these notes are as follows. We begin in £[2] with a brief thermal history of the CMB. 
We discuss the temperature and polarization anisotropy and acoustic peaks from recombination in 
Sj3]|4] and conclude in ^5} We take units throughout with % = c = ks = 1 and illustrate effects 
in the standard cosmological constant cold dark matter universe with adiabatic inflationary initial 
conditions (ACDM). 



2 Brief Thermal History 



In this section, we discuss the major events in the thermal history of the CMB. We begin in {2.1 



with the formation of the light elements and the original prediction of relic radiation. We continue 
in {2.2 with the processes that thermalize the CMB into a blackbody. Finally in {2.3 we discuss the 
recombination epoch where the main sources of temperature and polarization anisotropy lie. 



2.1 Nucleosynthesis and Prediction of the CMB 

Let us begin our brief thermal history with the relationship between the CMB and the abundance 
of light elements established at an energy scale of 10 2 keV, time scale of a few minutes, temperature 

1 



2 




of 10 9 K and redshift of z ~ 10 8 — 10 9 . This is the epoch of nucleosynthesis, the formation of the 
light elements. The qualitative features of nucleosynthesis are set by the low baryon-photon number 
density of our universe. Historically, the sensitivity to this ratio was used by Gamow and collaborators 
in the late 1940's to predict the existence and estimate the temperature of the CMB. Its modern use 
is the opposite: with the photon density well measured from the CMB spectrum, the abundance of 
light elements determines the baryon density. 

At the high temperature and densities of nucleosynthesis, radiation is rapidly thermalized to a 
perfect black body and the photon number density is a fixed function of the temperature n 7 oc T 3 
(see below). Apart from epochs in which energy from particle annihilation or other processes is 
dumped into the radiation, the baryon-photon number density ratio remains constant. 

Likewise nuclear statistical equilibrium, while satisfied, makes the abundance of the light elements 
of mass number A follow the expectations of a Maxwell-Boltzmann distribution for the phase space 
occupation number 

f A = e -( m A~VA)/T e -p 2 A /2m A T ^ /j\ 

where pa is the particle momentum, tua is the rest mass, and ha is the chemical potential. Namely 
their number density 



d 3 PA 
m A T, 



= gA{^f 2 e^-^l T . (2) 

Here gA is the degeneracy factor. In equilibrium, the chemical potentials of the various elements are 
related to those of the proton and neutron by 

fj, A = Z(x P + {A - Z)nn , (3) 



CMB Theory from Nucleosynthesis to Recombination 



3 



where Z is the charge or number of protons. Using this relation, the abundance fraction 

A-l 



X A = A^ = A^g A T 



A 



n b 



2vrT\ 3 / 2 2C(3)r ?fe7 



m b J 



B A /T X Z X A-Z 



(4) 



where X p and X n are the proton and neutron abundance, £(3) 
density. The two controlling quantities are the binding energy 



1.202, and n b is the baryon number 



Ba = Zm p + (A — Z)m n - tua , 



(5) 



and the baryon-photon number density ratio 7#, 7 = n b /n^. That the latter number is of order 10~ 9 
in our universe means that light elements form only well after the temperature has dropped below 
the binding energy of each species. Nuclear statistical equilibrium holds until the reaction rates drop 
below the expansion rate. At this point, the abundance freezes out and remains constant. 



Gamow's back of the envelope estimate (Gamow 1948 refined by Alpher and Herman 19481 was 



to consider the neutron capture reaction that forms deuterium 

n + p <-> d + 7 
with a binding energy of B 2 = 2.2MeV. Given Eqn. Q 



X 2 = A ( A ^'\ bi a?,)e B -l T X p X n 
it 2 \m b J 



(6) 



(7) 



a low baryon-photon ratio, and X n ~ X p ~ 1/2 for estimation purposes, the critical temperature 
for deuterium formation is T ~ 10 9 K. In other words the low baryon-photon ratio means that there 
are sufficient numbers of photons to dissociate deuterium until well below T ~ B 2 . Note that this 
condition is only logarithmically sensitive to the exact value of the baryon-photon ratio chosen for 
the estimate and so the reasoning is not circular. 

Furthermore, that we observe deuterium at all and not all helium and heavier elements means 
that the reaction must have frozen out at near this temperature. Given the thermally averaged cross 
section of 

(av) m 4.6 x KT^cmV 1 , (8) 

the freezeout condition 



n b (av) « H t 



(9) 



and the time-temperature relation i(T = 10 9 K) ~ 180s from the radiation dominated Friedmann 
equation, we obtain an estimate of the baryon number density 

n b (T = 10 9 K) 



1.2 x 10 17 cm -3 , 



(10) 



Comparing this density to the an observed current baryon density and requiring that n b oc a~ 3 and 
T oc a -1 yields the current temperature of the thermal background. For example, taking the modern 
0.02 and n b (a = 1) = 2.2 x 10~ 7 cm~ 3 yields the rough estimate 



value of £l b h 2 



2.2 x 10~ 7 cnr 
T(a = 1) r 



12K. 



:id 



This value is of the same order of magnitude as the observed CMB temperature 2.725K as well as 



the original estimates (Alpher and Herman [1948 Dicke et al. 19651. Modern day estimates of the 
baryon-photon ratio also rely on deuterium (e.g. Tytler et al.||2000 ). We shall see that the CMB has 
its own internal measure of this ratio from the acoustic peaks. Agreement between the nucleosynthesis 
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Fig. 2. Comptonization process. Energy injected into the CMB through heating of the electrons is thermalized 
by Compton scattering. Under the Kompaneets equation, a y-distortion first forms as low frequency photons 
gain energy from the electrons. Aft er multiple scatterings the distribution is thermalized to a chemical potential 
or /x-distortion. Adapted from Hu (1995). 



and acoustic peak measurements argue that the baryon-photon ratio has not changed appreciably 
since z ~ 10 s . 



2.2 Thermalization and Spectral Distortions 

Between nucleosynthesis and recombination, processes that create and destroy photons and hence 
thermalize the CMB fall out of equilibrium. The lack of spectral distortions in the CMB thus 
constrains any process that injects energy or photons into the plasma after this epoch. 

In a low baryon-photon ratio universe, the main thermalization process is double, also known as 
radiative, Compton scattering 

e~+7^e~+7 + 7. (12) 

The radiative Compton scattering rate becomes insufficient to maintain a blackbody at a redshift of 
( Danese and de Zotti|[l982 ) 



z therm = 2.0 x 10 6 (1 - r p /2)- 2 / 5 {^\ ' (13) 

corresponding to a time scale of order a few months. 

After this redshift, energy or photon injection appears as a spectral distortion in the spectrum of 
the CMB. The form of the distortion is determined by Compton scattering 

e~ + 7 <-> e~ + 7 , (14) 



since it is still sufficiently rapid compared with respect to the expansion while hydrogen remains 
ionized. Because a blackbody has a definite number density of photons at a given temperature, 
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Fig. 3. Low frequency thermalization by bremsstrahlung. Creation and absorption of photons at low frequencies 
by bremsstrahlung bring a chemical potential distortion back to blackbody at the electron temperature T e 
between zk and recombination. Adapted from Hu and Silk (19931. 



energy exchange via Compton scattering alone can only produce a Bose-Einstein spectrum for the 
photon distribution function 

f = e {E-v)/T _ I ( 15 ) 

with fi determined by the conserved number density and temperature. The evolution to a Bose- 



Einstein distribution is determined by solving the Kompaneets equation (Zeldovich and Sunyaev 



1969|. The Kompaneets equation is Boltzmann equation for Compton scattering in a homogeneous 
medium. It includes the effects of electron recoil and the second order Doppler shift that exchange 
energy between the photons and electrons. If energy is injected into the plasma to heat the electrons, 
Comptonization will try to redistribute the energy to the CMB. Consequently the first step in Comp- 
tonization is a so called "y distortion" in the spectrum as low energy photons in the Rayleigh Jeans 
tail of the CMB gain energy from the second order Doppler shift (see Fig. [2J. 

After multiple scatterings, the spectrum settles into the Bose-Einstein form with a "/i distor- 
tion." The transition occurs when the energy-transfer weighted optical depth approaches unity, i.e. 
t k(zk) =4 / dtn e aTT e /m e = 1 where n e is the electron density ot is the Thomson cross section, T e 
is the electron temperature and m e is the electron mass 



zk 



5.1 x 10 4 (1 - Y p /2) 



-1/2 



'n b h 2 ' 

0.02 



-1/2 



(16) 



Above this redshift energy injection creates a fi distortion, below a y distortion. 
At very low frequencies, bremsstrahlung 

e^+p^e~+p + j (17) 

is still efficient in creating and absorbing photons. In Fig. [3] we show the how a /i-distortion continues 



G 



GHz 




frequency (cm 1 ) 



Fig. 4. CMB spectrum from FIRAS. No spectral distortions from a blackbody have been discovered to date. 



to evolve until recombination which brings the low frequency spectrum back to a blackbody but now 
at the electron temperature. 

The best limits to date are from COBE FIRAS from intermediate to high frequencies: |yu| < 9 x 10 -5 
and \y\ < 1.5 x 10~ 5 at 95% confidence (see Fig [4] and Fixsen et al. 1996). After subtracting out 
galactic emission, no spectral distortions of any kind are detected and the spectrum appears to be a 
perfect blackbody of f = 2.725 ± 0.002K dMather et al.||l999|). 



2.3 Recombination 



While the recombination process 



p + e" 



H + 1 . 



(18) 



is rapid compared to the expansion, the ionization fraction obeys an equilibrium distribution just 
like that considered for light elements for nucleosynthesis or the CMB spectrum thermalization. As 
in the former two processes, the qualitative behavior of recombination is determined by the low 
baryon-photon ratio of the universe. 

Taking number densities of the Maxwell-Boltzmann form of Eqn. ([2]), we obtain 

Itt J 



-B/T 



n H 



where B 



m p + m e — ran = 13.6eV is the binding energy and we have set g p = g e = \gn 
the vanishingly small chemical potential of the photons, \i v + [i e = in equilibrium. 
Next, defining the ionization fraction for a hydrogen only plasma 



2. Given 



n, 



X e flf) . 



rib - n P = (1 - x e )n b , 



(20) 
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Fig. 5. Hydrogen recombination in Saha equilibrium vs. the calibrated 2-level calculation of RECFAST. The 
features before hydrogen recombination are due to helium recombination. 



we can rewrite Eqn. ( 19 ) as the Saha equation 

n en P = xj = 1 /m e T\ 3 / 2 ^ B/T 
nHn b 1 - x e n h \ 2-n ) 

Recombination occurs at a temperature substantially lower than T = B again because of the low 
baryon-photon ratio of the universe. We can see this by rewriting the Saha equation in terms of the 
photon number density 

x l = -B/kT * 1/2 ( m e c 2 \ 3/2 
1 - x e 25/2 %7C ( 3 ) \kT ) ■ { } 

Because rjby ~ 10™ 9 , the Saha equation implies that the medium only becomes substantially neutral 
at a temperature of T ~ 0.3eV or at a redshift of ~ 10 3 . At this point, there are not enough 
photons in the even in the Wien tail above the binding energy to ionize hydrogen. We plot the Saha 
solution in Fig. [5] 

Near the epoch of recombination, the recombination rates become insufficient to maintain ionization 
equilibrium. There is also a small contribution from helium recombination that must be added. The 
current standard for following the non-equilibrium ionization history is RECFAST ( Seager et al.|2000 ) 
which employs the traditional two-level atom calculation of [Peebles (19681 but alters the hydrogen 
case B recombination rate cub to fit the results of a multilevel atom. More specifically, RECFAST 
solves a coupled system of equations for the ionization fraction x% in singly ionized hydrogen and 
helium (i = H, He) 

dxj _ aBCjUHp f , _ x . ) _ x . x] (23) 
a ma H 

where njjp = (1 — Y p )rib is the total hydrogen plus proton number density accounting for the helium 



8 



mass fraction Y p , x e = n e jnu v = J2 x i 1S the total ionization fraction, n e is the free electron density, 
^max is the maximum X{ achieved through full ionization, 



< 

n>Hp 



-B ls /n 



c- 1 



1 + 



5rat 



Pa B e~ B2a/Tb 

A a + A2sls 

V 2?r 



(24) 



with g ra t the ratio of statistical weights, T& the baryon temperature, Bl the binding energy of the Lth 
level, A a the rate of redshifting out of the Lyman-o line corrected for the energy difference between 
the 2s and 2p states 



A-r 



7T 



and A 2s i s as the rate for the 2 photon 2s - 
B 2 s = B 2p = B ls /4, A 2sls = 8.22458s- 1 , 



(^max XijlT'Hp 

Is transition. For reference, for hydrogen B\ 
g Iat = 1, x max = 1. For helium B ls 



(25) 



3.967eV, B 2p = 3.366eV, A 2s i s = 51.3s" 1 , <? rat = 4. 



max 
3?max = 



13.598eV, 
24.583eV, B 2s = 



Y P /[A(1-Y P )]. 



Note that if the recombination rate is faster than the expansion rate asCjn^ p / H S> 1, the ion- 
ization solutions for Xi reach the Saha equilibrium dxi/dlna = 0. In this case s(x max — Xi) = XiX, 
or 

1 



\J (%ei ~\~ s) 2 ~\~ 4sx max (x e i -\- s) 



(26) 



1 



~\~ 3-max 
S 



1 



+ 2x max 
s 



+ 



where x e i = x e — Xi is the ionization fraction excluding the species. The recombination of hydrogenic 
doubly ionized helium is here handled purely through the Saha equation with a binding energy of 



1/4 the Bis of hydrogen and x T 



Y p /[4(1 — Yp)]. The case B recombination coefficients as 



a function of are given in Seager et al. (20001 as is the strong thermal coupling between T b and 
Tcmb- The multilevel-atom fudge that RECFAST introduces is to replace the hydrogen «b ^ 1.14as 
independently of cosmology. While this fudge suffices for current observations, which approach the 
~ 1% level, the recombination standard will require improvement if CMB anisotropy predictions are 
to reach an accuracy of 0.1% (e.g. |Switzer and Hirata|[2007l [Wong et al.||2007| ). 

The phenomenology of CMB temperature and polarization anisotropy is primarily governed by the 
redshift of recombination z* = a~ l — 1 when most of the contributions originate. This redshift though 
carries little dependence on standard cosmological parameters. This insensitivity follows from the 
fact that recombination proceeds rapidly once Bi s /Tb has reached a certain threshold as the Saha 
equation illustrates. Defining the redshift of recombination as the epoch at which the Thomson 
optical depth during recombination (i.e. excluding reionization) reaches unity, r rec (a*) = 1, a fit to 
the recombination calculation gives p!ul[2005] 

0.0105 / , 9 \ -0.028 

1 , (27, 



1089 



0.14 



0.024 



around a fiducial model of Q m h 2 = 0.14 and Q^h 2 = 0.024. 

The universe is known to be reionized at low redshifts due to the lack of a Gunn-Peterson trough in 
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Fig. 6. From temperature maps to power spectrum. The original temperature fluctuation map (top left) 
corresponding to a simulation of the power spectrum (top right) can be band filtered to illustrate the power 
spectrum in three characteristic regimes: the large-scale gravitational regime of COBE, the first acoustic peak 
where most o f the power lies, and the damping tail where fluctuations are dissipated. Adapted from Hu and 
White! pOQi). 



quasar absorption spectra. Moreover, large angle CMB polarization detections (see Fig. 19 ) suggest 
that this transition back to full ionization occurred around z ~ 10 leaving an extended neutral period 
between recombination and reionization. 



3 Temperature Anisotropy from Recombination 

Spatial variations in the CMB temperature at recombination are seen as temperature anisotropy by 
the observer today. The temperature anisotropy of the CMB was first detected in 1992 by the COBE 



10" 



across 



DMR instrument ( Smoot et al.|l99"2 ). These corresponded to variations of order AT/T 
10° - 90° on the sky (see Fig. [§. 

Most of the structure in the temperature anisotropy however is associated with acoustic oscillations 
of the photon-baryon plasma on ~ 1° scales. Throughout the 1990's constraints on the location of the 
first peak steadily improved culminating with the determinations of the TOCO ( |Miller et al. 1999), 
Boomerang, (de Bernardis et al. 2000) and Maxima-1 (Hanany et al. 2000) experiments. Currently 
from WMAP (Spergel et al. 2007) and ground based experiments, we have precise measurements of 
the first five acoustic peaks (see Fig. [7]). Primary fluctuations beyond this scale (below ~ 10') are 



damped by Silk damping (Silk 1968) as verified observationally first by the CBI experiment (Padin 



et al. 2001). 



In this section, we deconstruct the basic physics behind these phenomena. We begin with the 
geometric projection of temperature inhomogeneities at recombination onto the sky of the observer 
in {3.1 We continue with the basic equations of fluid mechanics that govern acoustic phenomena in 
£3.2 We add gravitational (£j3.3[), baryonic (jj3.4[), matter-radiation (|3.5|, and dissipational (£3.6) 
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Fig. 7. Temperature power spectrum from recent measurements from WMAP and ACBAR along with the best 
fit ACDM model. The main features of the temperature power spectrum includi ng the first 5 acoustic peaks 
and damping tail have now been measured. Adapted from Reichardt et al. (2008). 



effects in the sections that follow. Finally we put these pieces back together to discuss the information 
content of the acoustic peaks in §3.7| 



3. 1 Anisotropy from Inhomogeneity 

Given that the CMB radiation is blackbody to experimental accuracy (see Fig.Q, one can characterize 
its spatial and angular distribution by its temperature at the position x of the observer in the direction 
n on the observer's sky 

/(^,n,x) = [exp(2vr^/T(n;x) - l]" 1 , (28) 

where v = E/2ir is the observation frequency. The hypothetical observer could be an electron in the 
intergalactic medium or the true observer on earth or L2. When the latter is implicitly meant, we will 
take x = 0. We will occasionally suppress the coordinate x when this position is to be understood. 

For statistically isotropic, Gaussian random temperature fluctuations a harmonic description is 
more efficient than a real space description. For the angular structure at the position of the observer, 
the appropriate harmonics are the spherical harmonics. These are the eigenfunctions of the Laplace 
operator on the sphere and form a complete basis for scalar functions on the sky 

T(n) - f 



9(n) 



T 



Cm 



@£m^m(n) . 



(29) 



For statistically isotropic fluctuations, the ensemble average of the temperature fluctuations are 
described by the power spectrum 



(30) 



Moreover, the power spectrum contains all of the statistical information in the field if the fluctuations 
are Gaussian. Here Ci is dimensionless but is often shown with units of squared temperature, e.g. 
/iK 2 , by multiplying through by the background temperature today T . The correspondence between 
angular size and amplitude of fluctuations and the power spectrum is shown in Fig. [6} 

Let us begin with the simple approximation that the temperature field at recombination is isotropic 
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Fig. 8. From inhomogeneity to anisotropy. Temperature inhomogeneities at recombination are viewed at a 
distance D* as anisotropy on the observer's sky. 



but inhomogeneous and the anisotropy viewed at the present is due to the observer seeing different 
portions of the recombination surface in different directions (see Fig. pp. We will develop this picture 
further in the following sections with gravitational redshifts, dipole or Doppler anisotropic sources, 
the finite duration of recombination, and polarization. 

Under this simple instantaneous recombination approximation, the angular temperature fluctuation 
distribution is simply a projection of the spatial temperature fluctuation 

8(n) = j dDe(x)6(D - D*) , (31) 

where D = J dz/H is the comoving distance and denotes the distance a CMB photon travels 
from recombination. Here O(x) = (T(x) — T)/T is the spatial temperature fluctuation at recombi- 
nation. Note that the cosmological redshift does not appear in the temperature fluctuation since the 
background and fluctuation redshift alike. 

The spatial power spectrum at recombination can be likewise be described by its harmonic modes. 
In a flat geometry, these are Fourier modes 

d 3 k 



e(x)= y (2^ e(k)etk ' x - (32) 

More generally they are the eigenfunctions of the Laplace operator on the three dimensional space 



with constant curvature of a general Friedmann-Robertson- Walker metric (see e.g. Hu 2003). For 
a statistically homogeneous spatial distribution, the two point function is described by the power 
spectrum 

(0(k)*e(k')> = (2vr) 3 5(k - k')P(k) . (33) 

These relations tell us that the amplitude of the angular and spatial power spectra are related. 
A useful way of establishing this relation and quantifying the power in general is to describe the 
contribution per logarithmic interval to the variance of the respective fields: 

<e(x)9(x)> = J ^P(fc) = J dlnfc^P = J dlnkA 2 T (k). (34) 



12 

A scale invariant spectrum has equal contribution to the variance per e-fold A T = k 3 P(k) /2ir 2 = 
const. To relate this to the amplitude of the angular power spectrum, we expand equation (31) in 
Fourier modes 

©(n) = / ^3 ©(ky k ^ fi ■ (35) 
The Fourier modes themselves can be expanded in spherical harmonics with the relation 



4vr ]T i e j e (kD*)Y e * m (k)Y em (n) , (36) 

tm 



where jg is the spherical Bessel function. Extracting the multipole moments, we obtain 

Qim = J ^e(k)4Tri e j e (kD*)Y im (k) . (37) 



We can then relate the angular and spatial two point functions ( 30 ) and ( 33 ) 

(<d* em <di> m /) = 5 u >5 mm >4ir J din k j}{kD*)A%(k) = 8wS mm 'Ci. (38) 

Given a slowly varying, nearly scale invariant spatial power spectrum we can take A T out of the 
integral and evaluate it at the peak of the Bessel function kD* ~ I. The remaining integral can be 
evaluated in closed form J °° j\{x)d In x = l/[2£(£ + 1)] yielding the final result 

2ir 

Ce *W^i) A TW D *)- (39) 



Likewise even slowly varying features like the acoustic peaks of §3.2| also mainly map to multipoles 
of £ ~ kD* due to the delta function like behavior of jf. 

It is therefore common to plot the angular power spectrum as 

fc=^±MQ«A?.. (40) 

1/2 

It is also common to plot Cg ~ At, the logarithmic contribution to the rms of the field, in units of 
(jX. 

Now let us compare this expression to the variance per log interval in multipole space: 



(T(n)T(n)) = £ £ (T^T^Y^Y^,^) 



im I'm,' 

of J- 1 

E o E YL(m m (n) = e . (4i) 

em t 



For variance contributions from i 3> 1, 



E ^T°i « / ^ ^ ^ « / dln^^— . (42) 



47T J 4-7T J 2-K 

Thus is also approximately the variance per log interval in angular space as well. 
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3.2 Acoustic Oscillation Basics 

Thomson Tight Coupling — To understand the angular pattern of temperature fluctuations seen by 
the observer today, we must understand the spatial temperature pattern at recombination. That in 
turn requires an understanding of the dominant physical processes in the plasma before recombination. 

Thomson scattering of photons off of free electrons is the most important process given its relatively 
large cross section (averaged over polarization states) 

G T = -r— ^ = 6.65 x l(T 25 cm 2 . (43) 
3m e 

The important quantity to consider is the mean free path of a photon given Thomson scattering 
and a medium with a free electron density n e . Before recombination when the ionization fraction 
x e ~ 1 this density is given by 

n e = (1 - Y p )x e n b 

= 1.12 x 1(T 5 (1 - Y p )x e n b h 2 (l + z) 3 cm~ 3 . (44) 

The comoving mean free path Ac is given by 

A^ 1 = f = n e axa , (45) 

where the extra factor of a comes from converting physical to comoving coordinates and Y p is the 
primordial helium mass fraction. We have also represented this mean free path in terms of an 
scattering absorption coefficient f where dots are conformal time rj = / alt /a derivatives and r is the 
optical depth. 

Near recombination [z ~ 10 3 , x e ~ 1) and given n b h 2 ps 0.02, and Y p = 0.24, the mean free path is 

A c = - ~ 2.5Mpc . (46) 
r 

This scale is almost two orders of magnitude smaller than the horizon at recombination. On scales 
A 3> Ac photons are tightly coupled to the electrons by Thomson scattering which in turn are tightly 
coupled to the baryons by Coulomb interactions. 

As a consequence, any bulk motion of the photons must be shared by the baryons. In fluid 
language, the two species have a single bulk velocity t> 7 = v b and hence no entropy generation or 
heat conduction occurs. Furthermore, the shear viscosity of the fluid is negligible. Shear viscosity is 
related to anisotropy in the radiative pressure or stress and rapid scattering isotropizes the photon 



distribution. This is also the reason why in equation (31) we took the photon distribution to be 
isotropic but inhomogeneous. 

We shall see that the fluid motion corrects this by allowing dipole I = 1 anisotropy in the dis- 
tribution but no higher t modes. It is only on scales smaller than the diffusion scale that radiative 
viscosity (£ = 2) and heat conduction becomes sufficiently large to dissipates the bulk motions of the 
plasma (see £|3.6[). 



Zeroth Order Approximation — To understand the basic physical picture, let us begin our discussion 
of acoustic oscillations in the tight coupling regime with a simplified system which we will refine as 
we progress. 

First let us ignore the dynamical impact of the baryons on the fluid motion. Given that the baryon 
and photon velocity are equal and the momentum density of a relativistic fluid is given by (p + p)v, 
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where p is the pressure, and v is the fluid velocity, this approximation relates to the quantity 

R _ {Pb + Pb)vb = Pb+Pb = 3p b 

(47) 

where we have used the fact that p 7 oc T 4 so that its value is fixed by the redshifting background 
T = 2.725(l+z)K. Neglect of the baryon inertia and momentum only fails right around recombination. 

Next, we shall assume that the background expansion is matter dominated to relate time and scale 
factor. The validity of this approximation depends on the matter- radiation ratio 

Pr \ 0.15 ) VlO-V ' V ' 

and is approximately valid during recombination and afterwords. One expects from these arguments 
that order unit differences between the real universe and our basic description will occur. We will in 
fact use these differences in the following sections to show how the baryon and matter densities are 
measured from the acoustic peak morphology. 

Finally, we shall consider the effect of pressure forces and neglect gravitational forces. While this 
is not a valid approximation in and of itself, we shall see that for a photon-dominated system, the 
error in ignoring gravitational forces exactly cancels with that from ignoring gravitational redshifts 



that photons experience after recombination (see £3.3). 



Continuity Equation — Given that Thomson scattering neither creates nor destroys photons, the 
continuity equation implies that the photon number density only changes due to flows into and out 
of the volume. In a non expanding universe that would require 

n 7 + V • (n 7 v 7 ) = . (49) 

Since n 7 is the number density of photons per unit physical (not comoving) volume, this equation 
must be corrected for the expansion. The effect of the expansion can alternately be viewed as that 
of the Hubble flow diluting the number density everywhere in space. Because number densities scale 
as n 7 oc a~ 3 , the expansion alters the continuity equation as 

n 7 + 3n 7 — h V • (n 7 v 7 ) = . (50) 

Since we are interested in small fluctuations around the background, let us linearize the equations 
n 7 ~ n 7 + Srij and drop terms that are higher than first order in 5n 7 /n 7 and t> 7 . Note that t> 7 is first 
order in the number density fluctuations since as we shall see in the Euler equation discussion below 
it is generated from the pressure gradients associated with density fluctuations. 



The continuity equation (50) for the fluctuations becomes 

/ 5n. 



"7 



V-v 7 . (51) 



Since the number density n 7 oc T 3 , the fractional density fluctuation is related to the temperature 
fluctuation G as 

^ = 3f S 3 6. (52) 



CMB Theory from Nucleosynthesis to Recombination 



15 



Expressing the continuity equation in terms of G, we obtain 



e = -I V -v 7 . (53) 



Fourier transforming this equation, we get 

e = -^k-v 7 , (54) 
for the relationship between the Fourier mode amplitudes. 

Euler Equation — Now let us examine the origin of the fluid velocity. In the background, the velocity 
vanishes due to isotropy. However, Newtonian mechanics dictates that pressure forces will generate 
particle momentum as q = F. Newton's law must also be modified for the expansion. If we 
associate the de Broglie wavelength with the inverse momentum, this wavelength also stretches with 
the expansion. For photons, this accounts for the redshift factor. For non-relativistic matter, this 
means that bulk velocities decay with the expansion. In either case, we generalize Newton's law to 
read 

q+% = F. (55) 
For a collection of particles, the relevant quantity is the momentum density 

/d?q 

and likewise the force becomes a force density. For the photon-baryon fluid, this force density is 
provided by the pressure gradient. The result is the Euler equation 

[(p 7 +p 7 )v 7 ]" = -4^(p 7 +p 7 )v 7 - Vp 7 , (57) 

where the 4 on the right hand side comes from combining the redshifting of wavelengths with that of 
number densities. Since photons have an equation of state p 7 = p 7 /3, the Euler equation becomes 

4 • Ir, 

-p 7 v 7 = --Vp 7 , 

v 7 = -Ve. (58) 

In Fourier space, the Euler equation becomes 

v 7 = -ike . (59) 

The factor of i here represents the fact that the temperature maxima and minima are zeros of the 
velocity in real space due to the gradient relation, i.e. temperature and velocity have a n/2 phase 
shift. It is convenient therefore to define the velocity amplitude to absorb this factor 

v 7 = — iv-fk . (60) 

The direction of the fluid velocity is always parallel to the wavevector k for linear (scalar) perturba- 
tions and hence we can write the Euler equation as 

v 1 = kQ. (61) 
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Fig. 9. Acoustic oscillation basics. All modes start from the same initial epoch with time denoted by the 
sound horizon relative to the sound horizon at recombination s*. (a) Wavenumbers that reach extrema in 
their effective temperature + 9 (accounting for gravitational redshifts ^3.3 1 at s* form a harmonic series 
k n = mr/s*. ( b) Amplitude of the fluct uations is the same for the maxima and minima without baryon inertia. 
Adapted fromlHu and Dodelsonl (120021). 



Acoustic Peaks — Combining the continuity (54) and Euler (61) equations to eliminate the fluid 
velocity, we get the simple harmonic oscillator equation 

9 + c^ 2 9 = 0, (62) 

where the adiabatic sound speed c 2 . = 1/3 for the photon-dominated fluid and more generally is 
defined as 

^ ■ (63) 
Pi 

The solution to the oscillator equation can be specified given two initial conditions 0(0) and u 7 (0) 
or 8(0), 

8(77) = 9(0) cos(£;s) + sin(£;s) , (64) 

KC S 

where the sound horizon is denned as 

'a.*,. (65) 



In real space, these oscillations appear as standing waves for each Fourier mode. 

These standing waves continue to oscillate until recombination. At this point the free electron 
density drops drastically (see Fig. [5]) and the photons freely stream to the observer. The pattern of 
acoustic oscillations on the recombination surface seen by the observer becomes the acoustic peaks 
in the temperature anisotropy. 

Let us focus on the adiabatic mode which starts with a finite density or temperature fluctuation 
and vanishing velocity perturbation. At recombination 77*, the oscillation reaches (see Fig. [9]) 

9(77*) = 9(0) cos(fcs*). (66) 

Considering a spectrum of k modes, the critical feature of these oscillations are that they are tempo- 
rally coherent. The underlying assumption is that fluctuations of all wavelengths originated at 77 = 
or at least 77 <C 77*. Without inflation this would violate causality for long wavelength fluctuations, i.e. 
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the analogue of the horizon problem for perturbations. With inflation, super horizon modes originate 
during an inflationary epoch rfc <C ry*. 

Modes caught in the extrema of their oscillation follow a harmonic relation 

k n s* = mr , n = 1,2,3... (67) 

yielding a fundamental scale or frequency, related to the inverse sound horizon 

kA = 7r/s* . (68) 

Since the power spectrum is proportional to the square of the fluctuation, both maxima and minima 
contribute peaks in the spectrum. Observational verification of this harmonic series is the primary 
evidence for inflationary adiabatic initial conditions (Hu and White 1996). 

The fundamental physical scale is translated into a fundamental angular scale by simple projection 
according to the angular diameter distance Da 

e A = x A /D A , 

l A = k A D A , (69) 



(see Eqn. (39)). In a flat universe, the distance is simply Da = D = fyo — ?/* ~ VOi the horizon 
distance, and kA = vr/s* = \/37r/ry* so 

9a « ^ ■ (70) 
Vo 

Furthermore, in a matter-dominated universe rj tx a 1 / 2 so 6a ~ 1/30 ~ 2° or 

l A « 200 . (71) 



We shall see in £3.7 that radiation and dark energy introduce important corrections to this prediction 
from their influence on 77* and D* respectively. Nonetheless it is remarkable that this simple argument 
predicts the basics of acoustic oscillations: their existence, coherence and fundamental scale. 

Acoustic Troughs: Doppler Effect — Acoustic oscillations also imply that the plasma is moving 
relative to the observer. This bulk motion imprints a temperature anisotropy via the Doppler effect 



= n • v 7 . (72) 

dop 



Averaged over directions 



/ATN 



which given the acoustic solution of Eqn. (|64|) implies 



(73) 



% = -^Q = ^kc s @(0) S m(ks) 
V3 k k 

= 9(0)sin(fcs). (74) 

Interestingly, the Doppler effect for the photon-dominated system is of equal amplitude and tt/2 
out of phase: extrema of temperature are turning points of velocity. If we simply add the fc-space 
temperature and Doppler effects in quadrature we would obtain 

f^T) 2 = 9 2 (0)[cos 2 (A;s) +sin 2 (£;s)] = 9 2 (0) . (75) 
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Fig. 10. Doppler effect. The Dopplcr effect provides fluctuations of comparable strength to the local tempera- 
ture fluctuations from acoustic oscillations in fc-space (lower) providing features at the troughs of the latter. In 
angular space, projection effects smooth the Doppler features leaving an acoustic morph ology that reflects the 
temperature oscillations. The peak height modulatio n co mes from the baryon inertia (j |3.4| and the gradual 
increase in power with I from radiation domination (j|3.5[). 



In other words there are no preferred /c-modes at harmonics and a scale invariant initial temperature 
spectrum would lead to a scale invariant spatial power spectrum at recombination. However the 
Doppler effect carries an angular dependence that changes its projection on the sky n • v 7 oc n • k. In 
a coordinate system where z || k, this angular dependence yields an extra factor of Yio in the analogue 



of Eqn. (38). This extra factor can be reabsorbed into the total angular dependence through Clebsch- 



Gordan recoupling ( Hu and White| [l997) 



Y 1Q Y eo ^Y e±10 . (76) 

The recoupling implied for the radial harmonics changes ji{x) — * j'i(x). The projection kernel j[{x) 
lacks a strong feature at i ~ x and so Doppler contributions in k are spread out in £. This is 
simply a mathematical way of stating that the Doppler effect vanishes when the observer is looking 
perpendicular to v || k whereas it is in that direction that the acoustic peaks in temperature gain 
most of their contribution. The net effect, including baryonic effects that we discuss below, is that 
the peak structure is dominated by the local temperature at recombination and not the local fluid 
motion. 



3.3 Gravito- Acoustic Oscillations 

Thus far we have neglected gravitational forces and redshifts in our discussion of plasma motion. The 
true system exhibits gravito- acoustic or Jeans oscillations. We were able to employ this swindle to 
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get the basic properties of acoustic oscillations because in a photon-dominated plasma the effect of 
gravitational forces and gravitational redshifts exactly cancel in constant gravitational potentials. To 
go beyond the photon dominated plasma and matter dominated expansion approximation, we now 
need to include these effects. Furthermore, we shall see that the gravitational potential perturbations 
from inflation are also the source of the initial temperature fluctuation. 

Continuity Equation and Newtonian Curvature — The photon continuity equation for the number 
density is altered by gravity since the presence of a gravitational potential alters the coordinate 
volume. Formally in general relativity this comes from the space-space piece of the metric - a spatial 
curvature perturbation <£: 

ds 2 = o 2 [-(l + 2^)d7] 2 + (1 + 2$)dx 2 ] (77) 

for a flat cosmology. 

We can think of this curvature perturbation as changing the local scale factor a — * a(l + 3>) so that 
the expansion dilution is generalized to 

*-* + *. (78) 
a a 

Hence the full continuity equation is now given by 

(5n~Y = — 3<5n 7 3n 7 6 — n 7 V • v 7 , (79) 

a 

or 

S = - 1 kv 1 -<Z>. (80) 
3 

Euler Equation and Newtonian Forces — Likewise the gravitational force from gradients in the gravi- 
tational potential (^f, formally the time-time piece of the metric perturbation) modifies the momentum 
conservation equation. The Newtonian force F = -raV$ generalized to momentum density brings 
the Euler equation to 

u 7 = fc(0 + ¥). (81) 

General relativity says that $ and ^ are the relativistic analogues of the Newtonian potential and 
that $ ~ — ^ in the absence of sources of anisotropic stress or viscosity. 



Photon Dominated Oscillator — We can again combine the continuity equation (80) and Euler equa- 



tion (81 ) to form the forced simple harmonic oscillator system 

e + c 2£: 2 e = - — (82) 
3 

Note that the effect of baryon inertia is still absent in this system. To make further progress in un- 
derstanding the effect of gravity on acoustic oscillations we need to specify the gravitational potential 
v|/ « — $ from the Poisson equation and understand its time evolution. 

Poisson Equation and Constant Potentials — In our matter-dominated approximation, $ is generated 
by matter density fluctuations A m through the cosmological Poisson equation 

k 2 <5> = 4irGa 2 p m A m , (83) 



where the difference from the usual Poisson equation comes from the use of comoving coordinates 
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for k (a 2 factor), the removal of the background density into the background expansion (p m A m ) and 
finally a coordinate subtlety that enters into the definition of A m . In general, the relativistic Poisson 
equation would have contributions from the momentum density. Alternately we can relate 

A m = ^ + 3^ (84) 

to the density fluctuation in a coordinate system that comoves with the matter (formally through 
a gauge transformation). Beyond the matter dominated approximation, the Poisson equation would 
carry contributions from all of the species of energy density and the comoving coordinate system 
would also reflect the total. 

In a matter-dominated epoch (or in fact any epoch when gravitational potential gradients and not 
stress gradients dominate the momentum equation), the matter Euler equation implies v m ~ krj$?. 
The matter continuity equation then implies A m ~ —kr\v m ~ — {kr\) 2 ^. The Poisson equation then 
yields ~ A m /(kr]) 2 ~ — Here we have used the Friedmann equation H 2 = 8TiGp rn /3 and 
rj = J d\a.a/{aH) ~ \/{aH). In other words, the density perturbation A m grows at exactly the right 
rate to keep the gravitational potential constant. More formally, if stress perturbations are negligible 
compared with density perturbations ( bp <C bp ) then the potential will remain constant in periods 
where the background equation of state p/p is constant. With a varying equation of state, it is the 
comoving curvature C rather than the Newtonian curvature that is strictly constant ( Bardeen||1980 ). 



Effective Temperature — If the gravitational potential is constant, we can rewrite the oscillator 



equation (82) as 

@ + ^ + c 2 s k 2 (e + ^) = 0. (85) 



The solution in Eqn. ( 64 ) for adiabatic initial conditions then generalizes to 

[G + (77) = [9 + *] (0) cos{ks) . (86) 

Like a mass on spring in a constant gravitational field of the earth, the solution just represents 
oscillations around a displaced minimum. 

Furthermore, G + ^ is also the observed temperature fluctuation. Photons lose energy climbing 
out of gravitational potentials at recombination and so the observer at the present will see 

AT 

— = e + *. (87) 

Therefore from the perspective of the observer, the acoustic oscillations ignoring both gravitational 
forces and gravitational redshifts is unchanged: initial perturbations in the effective temperature 
oscillate around zero with a frequency given by the sound speed. What the consideration of gravity 
adds is a way of connecting the initial conditions to inflationary curvature fluctuations. 

The Sachs-Wolfe 1/3 — The effective temperature perturbation includes both the local temperature 
in Newtonian coordinates and the gravitational redshift factor ^f. For ks* <C 1, the oscillator is 



frozen at its initial conditions and the total is called the Sachs and Wolfe (1967) effect. The division 
into two pieces is actually an artifact of the coordinate system. Both pieces are determined by the 
initial curvature perturbation from inflation. 

To see this, let us relate both the Newtonian potential and local temperature fluctuation to a 
change in time coordinate from comoving coordinates to Newtonian coordinates (formally through 
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Fig. 11. Acoustic oscillations with gravitational forcing and dissipational damping. For a mode that enters 
the sound horizon during radiation domination, the gravitational potential decays at horizon crossing and 
drives the acoustic amplitude higher. As the photon diffusion length increases and becomes comparable to 
the wavelength, r adia tive viscosity Tr ~, is generated from qua drupole anisotropy leading to dissipation and 
polarization (see { 4.4 1. Adapted from Hu and Dodelson (20021. 



a gauge transformation, see White and Hu 1997 for a more detailed treatment). The Newtonian 
gravitational potential ^ is a perturbation to the temporal coordinate (see Eqn. (77)) 

St 

T 



Given the Friedmann equation, this is equivalent to a perturbation in the scale factor 



t 



da 
aH 



oc 



da 



ap 



1/2 



oc a 



3(l+w)/2 



(89) 



(90) 



where w = p/p. During matter domination w = and 

5a _ 2 St 
a 3 t 

Since the CMB temperature is cooling as T oc a -1 a local change in the scale factor changes the 
local temperature 

e = - (5a 

a 

Combining this with \& to form the effective temperature gives 



3 



(91) 



e + * = 



(92) 



The consequence is that overdense regions where ^ is negative (potential wells) are cold spots in the 
effective temperature. 

Inflation provides a source for initial curvature fluctuations. Specifically the comoving curvature 
perturbation £ becomes a Newtonian curvature of 



(93) 
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Fig. 12. Acoustic oscillations with baryons. Baryons add inertia to the photon-baryon plasma displacing the 
zero point of the oscillation and making compressional peaks (minima) larger than rarefaction peaks (maxima) . 
The absolute value of the fluctuation in effective temperature is shown in dotted lines. Adapted from Hu and 



Dodelson (2002). 



in the matter-dominated epoch. The initial amplitude of scalar curvature perturbations is usually 
given as As = 5^ which characterizes the variance contribution per e-fold to the curvature near some 
fiducial wavenumber k n (see Eqn. 34) 



A 2 = = a? 

C 2vr 2 C 



71- 1 



where n = 1 for a scale invariant spectrum. Combining these relations for n = 1 

i(£+l)Ci 



2tt 
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(94) 



(95) 



in the Sachs- Wolfe limit. The 10 5 fluctuations measured by COBE then correspond to 6^ ~ 5 x 10 5 . 



3.4 Baryonic Effects 

The next level of detail that we need to add is the inertial effect of baryons in the plasma. Since 



equation (47) says that the baryon momentum becomes comparable to the photon momentum near 
recombination, we can expect order unity corrections on the basic acoustic oscillation picture from 
baryons. With the precise measurements of the first and second peak, this change in the acoustic 
morphology has already provided the most sensitive measure of the baryon-photon ratio to date 
exceeding that of big bang nucleosynthesis. 

Baryon Loading — Baryons add extra mass to the photon-baryon plasma or equivalently an enhance- 
ment of the momentum density of the plasma given by R = (pb + Pb)/(Pj + P^)- Specifically the 
momentum density of the joint system 

(p 7 + p 7 )t; 7 + (p b + pb)v h = (1 + R)(p-y + p 7 )w 7 b (96) 

is conserved. For generality, we have introduced the baryon velocity Vb and the momentum-weighted 



velocity v^b but in the tightly coupled plasma Vb ~ v^b ~ w 7 (c./. ^3.6). 
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The Euler equation (81) becomes 

[(1 + R)(Pj +Pi)v 1 b\ = -4^(1 + -R)(p 7 +p 7 )v 7b 

-Vp 7 - (l + i?)(p 7 +.p 7 )V#. (97) 

This equation takes the same form as the photon-dominated system except for the (1 + R) terms 
which multiply everything but the pressure gradient terms since the pressure comes predominantly 
from the photons. We can rewrite the equation more compactly as 

[(I + R)v lb \ = k® + (I + R)m . (98) 



Oscillator with Gravity and Baryons — The photon continuity equation (80) remains the same so 
that the oscillator equation becomes 

[(1 + R)Q\ + \k 2 @ = -\k 2 {l + R)W - [(1 + R)$y . (99) 

This equation is the final oscillator equation for the tight coupling or perfect fluid regime. 

We can make several simplifications to illuminate the impact of baryons. First let us continue to 
use the matter dominated approximation where = — <I> =const. Next let us make the adiabatic 
approximation where the change in R is slow compared with the frequency of oscillation R/R <C u) = 
kc s . In that case, Eqn. (99) looks like an oscillator equation with a fractional change in the mass 
given by R and a change in the sound speed 

(100) 



2 _ Vl +Pb _ 1 



s ^ + p b 3(1 + R) ' 

Consequently, the solution in Eqn. ( |86[ ) is modified as 

[9 + (l + #)*](??) = [9 + (1 + R)^}{0) cos ks. (101) 

This solution is reminiscent of that of adding mass to the spring in a constant gravitational field of 
the earth. 

There are 3 effects of baryon loading. First the amplitude of oscillations increases by a factor of 
1 + 3R 

[0 + (l + J2)¥](O) = J(l + 3#)*(0). (102) 

Next the equilibrium point of the oscillation is now shifted so that relative to zero effective tempera- 
ture, the even and odd peaks have different amplitudes 

[6 + *]„ = [±(l + 3fl)-3-R]i*(0), 

[e + tf]i-[e + ¥] 2 = [-6J2]i*(0). (103) 



In particular, baryon loading increases the heights of the odd peaks over the even peaks (see Fig. 12 ). 
Finally, the lowering of the sound speed changes the acoustic scale in Eqn. (69) as 

I A oc Vl + R ■ (104) 

The effects of baryon loading in a full calculation are actually smaller since R is growing in time. 

Baryon-Photon Momentum Ratio Evolution — One can get a handle on the effect of evolution of R 
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by again equating the system to the analogous physical oscillator. The baryons add inertia or mass 
to the system and for a slowly varying mass the oscillator equation has an adiabatic invariant 

E 1 1 

— = -m eff uM 2 = -(1 + R)kc s A 2 oc A 2 (l + R) 1/2 = const. (105) 
u 2 2 

Amplitude of oscillation A oc (1 + -R) -1 / 4 decays adiabatically as the photon-baryon ratio changes. 
This offsets the gain in the overall amplitude from Eqn. (102). Coupled with uncertainties in the 
distance to recombination in interpreting the I a measurement, this leaves the modulation of the peak 
heights as the effect that provides most of the information about the baryon-photon ratio in the 
acoustic peaks. 



3.5 Matter- Radiation Ratio 

Next we want to go beyond the matter-dominated expansion approximation. The universe only be- 



comes matter dominated in the few e-folds before recombination (see Eqn. (48)). Peaks corresponding 
to wavenumbers that began oscillating earlier carry the effects of the prior epoch of radiation dom- 
ination. These effects come in through the evolution of the gravitational potential which acts as a 



forcing function on the oscillator through Eqn. (99). 



Potential Decay — The argument given in §3.3| for the constancy of the gravitational potential depends 
crucially on gravity being the dominant force affecting the total density. When radiation dominates 
the total density, radiation stresses become more important than gravity on scales smaller than the 
sound horizon. The total density fluctuation stops growing and instead oscillates with the acoustic 
frequency. The Poisson equation in the radiation dominated epoch 

k 2 $ = 4irGa 2 p r A r (106) 



then implies that <E> oscillates and decays with an amplitude oc a~ 2 (see Fig. 11). As an aside, the 
relativistic stresses of dark energy make the gravitational potential decay again during the acceleration 
epoch and lead to the so called integrated Sachs- Wolfe effect. 



Radiation Driving — An examination of Fig. 11 shows that the time evolution of the gravitational 



potential is in phase with the acoustic oscillations themselves and act a driving force on the acoustic 
oscillations. We can estimate the effect on the amplitude of oscillations in the limit that the force is 



fully coherent. In that case we can take the continuity equation (80) and simply integrate it 



[9 + *](?/) = [9 + *](0) +A^- A$ 

= ±*(0) - 2$(0) = §*(0) . (107) 

This estimate gives an acoustic amplitude that is 5x that of the Sachs- Wolfe effect. This enhancement 
only occurs for modes that begin oscillating during the radiation dominated epoch, i.e. the higher 
peaks. The net effect is a gradual ramp up of the acoustic oscillation amplitude across the horizon 



wavenumber at matter radiation equality (see Fig. 10). 



In fact the coherent approximation is exact for a photon-baryon fluid but must be corrected for 
the neutrino contribution to the radiation density. 

External Potential Approach — For pedagogical purposes it is sometimes useful to go beyond the 
coherent approximation. Neutrino corrections are one example; isocurvature initial conditions are 
another. 
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We have seen that the solutions to homogeneous equation for the oscillator equation ( 82 ) are 

{1 + R)- 1/4 cos{ks) , (1 + R)- 1/4 sm{ks) (108) 

in the adiabatic or high frequency limit. Considering the potentials as external, we can solve for the 
temperature perturbation as ( |Hu and Sugiyania 1995) 

(1 + R) 1 ' 4 ®^) = 9(0)cos(A;s) + ^ fe(0) + \R(0)Q(0) 



sin ks 



+ T Jo */(l + R') 3/4 sm[ks - ks']F(r)') , 



where 



R 



(109) 



;no) 



l + R 3 

By including the neutrino effects in the gravitational potential, we can show from this approach that 
radiation driving actually creates an acoustic amplitude that is close to 4x the Sachs- Wolfe effect. 



3.6 Damping 

The final piece in the acoustic oscillation puzzle is the damping of power beyond £ ~ 10 3 shown in 
Fig. [6} Up until this point, we have considered the oscillations in the tight coupling approximation 
where the photons and baryons respond to pressure and gravity as a single perfect fluid. Fluid 



imperfections are associated with the Compton mean free path in Eqn. (46). Dissipation becomes 



strong at the diffusion scale, the distance a photon can random walk in a given time rj (Silk 1968) 



9 = --t; 7 -i>, 6 b = -kv b - 3$ . (112) 



X D = VNXc = ^Jr]/X c Xc = Vv^c- (HI) 

This scale is the geometric mean between the horizon and mean free path. Given that Xo/r}* ~ few 
percent, we expect that the n > 3 peaks to be affected by dissipation. To improve on this estimate, 
we develop the microphysical description of dissipation next. 

Continuity Equations — To treat the photons and baryons as separate systems, we now need to 
supplement the photon continuity equation with the baryon continuity equation 

k 
3 

The baryon equation follows from number conservation with p b = m^rii) and 5 b = bp b j p b . 

Euler and Navier-Stokes Equations — The momentum conservation equations must also be separated 
into photon and baryon pieces 

k 

i) 7 = k(@ + <£) - -7r 7 - f(Vj - v b ) , 

v b = -^ v b + + f(u 7 - v b )/R, 
where the photons gain an anisotropic stress term 7r 7 from radiation viscosity. The baryon equation 



follows from the same derivation as in ^3.2 where the redshift of the momentum is carried by the 
bulk velocity instead of the redshifting temperature. Finally there is a momentum exchange term 
from Compton scattering. Note that the total momentum in the system is conserved and hence the 
scattering terms come with opposite sign. 



2G 

Viscosity — Radiative shear viscosity is equivalent to quadrupole moments in the temperature field. 
These quadrupole moments are generated by radiation streaming from hot to cold regions much like 
how temperature inhomogeneity are converted to anisotropy in Fig. |8j 

In the tight coupling limit where f/k, the optical depth through a wavelength of the fluctuation is 
high one therefore expects 



k 

V>y~ , 
T 



(113) 



since it must be generated by streaming and suppressed by scattering. A more detailed calculation 
from the Boltzmann or radiative transfer equation says flKaiser||l"983] ) 



f 



where A v = 16/15 once polarization effects are incorporated 

k k 

v = k(Q + $0 A) T D 7 . 

3 T 



(114) 



(115) 



The oscillator equation with viscosity becomes 



d 



2„2 



k z c 



<4-ir(cj 2 e) + ^A v e + k 2 cje 



k 



d 



' drj ' " t 

As in a mechanical oscillator, a term that depends on provides a dissipational term to the solutions. 

Heat Conduction — Relative motion between the photons and baryons also damp oscillations. By 
expanding the continuity and momentum conservation equations in the small number k jr one obtains 
for the full oscillator equation 



c 2 s ^(c; 2 Q) + ^[A V + A h }& + k z ci& 
dr/ t 



s dr, { s ' 



where 



A h 



R 2 



1 + R 



(116) 



Dispersion Relation — We can solve the damped oscillator equation in the adiabatic approximation 
by taking a trial solution O oc exp(i J ojdrj) to obtain the dispersion relation 

i kc.o 



1±- — (A V + A h ) 
z r 



;ii7 



The imaginary term in the dispersion relation gives an exponential damping of the oscillation ampli- 
tude 



exp(i / udrj) 



e ±iks exp[-k 2 I dri^(A v + A h 



lc 2 



2 f 



e ±iks e W [-(k/k D ) 2 } 



where the diffusion wavenumber is given by 



D 



dr]— 



1 



16 



+ 



R 2 



f 6(1 + 11) V 15 {1 + R) 
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Note that in both the high and low R limits 



Hence the dissipation scale is 



- < i2o) 

2TT 2tT .1.1/9 . 

AD = to ~ ^ > ' (121) 

and comparable to the geometric mean between horizon and mean free path as expected from the 
random walk argument . For a baryon density of Sl^/i 2 ~ 0.02, radiation viscosity is responsible for 
most of the dissipation and we show the correspondence between viscosity generation and dissipation 
in Fig. [TTJ 

Since the diffusion length changes rapidly through recombination and the medium changes from 
optically thick to optically thin, the damping estimates above are only qualitative. A full Boltzmann 
(radiative transfer solution) shows a more gradual, but still exponential, damping of roughly 

V t « expH^) 1 ' 25 ] , (122) 



with a damping scale of in = 2irD if /\r) and (Hu 20051 



Mpc ^ 0.14 J ^0.024 J v ' 

for small changes around the central values of Sl m /i 2 and VL^h?. 

This envelope also accounts for enhanced damping due to the finite duration of recombination. 



Instead of a delta function in the projection equation (31) we have the visibility function fe~ T that 
acts as a smearing out of any contributions with wavelengths shorter than the thickness of the 
recombination surface that survive dissipation. 



3. 7 Information from the Peaks 

In the preceding sections we have examined the physical processes involved in the formation of the 
acoustic peaks and explained their sensitivity to the energy content and expansion rate of the universe. 
Converting the measurements into parameter constraints of course requires a more accurate numerical 
description. Numerical codes that solve the Einstein- Boltzmann radiative transfer equations for the 
CMB and matter ( jPeebles and Yu|[l970{ |Bond and Efstathiou||1984| |Vittorio and Silk|[T984[ ) are now 
accurate at the ~ 1% level for the acoustic peaks for publically available codes (Seljak and Zaldarriaga 



1996 Lewis et al. 2000). Their numerical precision on the other hand is substantially better and 

10 3 . 



approaches the 0.1% level required for cosmic variance limited measurements out to 



The 



accuracy is now limited by the input physics, mainly recombination (see 2.3). In this section, we 
relate the qualitative discussion of the previous sections to the quantitative information content of 
the peaks. 

First Peak: Curvature and Dark Energy — The comparison between the predicted acoustic peak scale 
A,4 and its angular extent provides a measurement of the angular diameter distance to recombination. 
The angular diameter distance in turn depends on the spatial curvature and expansion history of the 
universe. 
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Fig. 13. Angular diameter distance and curvature. In a non-flat (here closed) universe, the apparent or angular 
diameter distance Da = Lda does not equal the radial distance traveled by the photon. Objects in a closed 
universe are further than they appear whereas in an open universe they are closer than they appear. 



Sensitivity to the expansion history during the acceleration epoch comes about through the radial 
distance a photon travels along the line of sight 

d - r — 

*~Jo H(z) 

With the matter and radiation energy densities measured, the remaining contributor to the expansion 
rate H (z) is the dark energy. 



(124) 



However in a curved universe, the apparent or angular diameter distance Da in Eqn. (69) is no 



longer the distance a photon travels radially along the line of sight. The radius of curvature of the 
space is given in terms of the total density fitot i n units of the critical density as 



R~ 



tot 



1) 



(125) 



A positive curvature space has O to t > 1 and a real radius of curvature. A negatively curved space 
has Otot < 1 an d an imaginary radius of curvature. The positively curved space is shown in Fig. 13 
The curvature makes a transverse distance L related to its angular extent da as L = daD^ with 



D A = Rsin(D/R) . 



(126) 



The same formula applies for negatively curved spaces but is more conveniently expressed with the 
relation 



Rsin(D/R) = \R\sinh(D/\R\ 



(127) 



for imaginary R. In a positively curved geometry Da < D and objects are further than they appear. 
In a negatively curved universe R is imaginary and Rsin(D/ R) = i\R\ s'm(D/i\R\) = \R\ sinh(D/|i?|) 
- and Da > D objects are closer than they appear. Since the detection of the first acoustic peak it 



has been clear that the universe is close to spatially flat (Miller et al. 1999 de Bernardis et al. 2000 



Hanany et al.pOOO ). How close and how well-measured D* is for dark energy studies depends on the 
calibration of the physical scale Xa = 2s*, i.e. the sound horizon at recombination. 
The sound horizon in turn depends on two things, 

2\/3 / 



In 



VI + R* + VR* + r *R* 

i + V r *^* 



(128) 
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Fig. 14. Curvature and dark energy. Given a fixed physical scale for the acoustic peaks (fixed Ut,h 2 and 
Q, m h 2 ) the observed angular position of the peaks provides a measure of the angular diameter distance and the 
parameters it depends on: curvature, dark energy density and dark energy equation of state. Changes at low 
I multipoles are due to the decay of the gravitational potential after matter domination from the integrated 
Sachs- Wolfe effect. 



the baryon-photon momentum density ratio 



R* 



3 pb 
4p 7 



0.729 



'n b h? s 
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a* 
10" 3 
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and the expansion rate prior to recombination which is determined by the matter radiation ratio 

-i 



Pr 

Pn 



0.297 




a* 



(130) 



The calibration of these two quantities involves the higher acoustic peaks. The bottom line is that 
the limiting factor in the calibration is the precision with which the matter density is known 

1 6(n m h 2 ) 

4 <rv} hi 



5 Da* 



Da* 



(131) 



In principle the Planck satellite can achieve a 1% measurement of the matter density or a ~ 0.25% 
measure of distance. Current errors on the distance are ~ 1-2%. 

Second Peak: Baryons — The baryon-photon ratio controls the even-odd modulation of peak heights 



through the baryon loading effect (see £3.4). The second peak represents rarefaction of the acoustic 
wave in a gravitational potential and hence is suppressed in amplitude by the baryon inertia. The 



dependence of the spectrum on the baryon density Q b h is shown in Fig. 15 Since the first tentative 



detection of the second peak (de Bernardis et al. 20001 the CMB has placed limits on the baryon 
density. Currently its measurement, largely from WMAP ( Spergel et al.|2007 ), gives n b h 2 = 0.0227± 
0.0006 (Reichardt et al. 20081. This is well enough constrained that associated errors on the sound 



horizon and shape of the matter power spectrum (see below) are small. 

Third Peak: Dark Matter — The third peak begins to show the effects of the matter-radiation ratio on 
the overall amplitude of the acoustic peaks. Furthermore, decay in the gravitational potential during 
radiation domination would reduce the baryon loading effect and change the peak height ratios of the 
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Fig. 15. Baryons and matter. Baryons change the relative heights of the even and odd peaks through their 
inertia in the plasma. The m atter-radiation rat i o also changes the overall amplitude of the oscillations from 
driving effects. Adapted from |Hu and Dodelson (2002). 



second and third peaks (e.g. Hu et aL|2001 ). The dependence of the spectrum on the baryon density 



O m /i 2 is shown in Fig. 15[ Constraints on the third peak from the DASI experiment (iPryke et al. 



2001 ) represented the first direct evidence for dark matter at the epoch of recombination. Current 



constraints from a combination of WMAP and higher resolution ground and balloon based data yield 
Qmh 2 = 0.135 db 0.007 (Reichardt et al. (2008)). Since this parameter controls the error on the 



distance to recombination through equation (131) and the matter power spectrum (see below), it is 



important to improve the precision of its measurement with the third higher peaks. 

Damping Tail: Consistency — Under the standard thermal history of ^2] and matter content, the 
parameters that control the first 3 peaks also determine the structure of the damping tail at I > 10 3 : 
namely, the angular diameter distance to recombination D*, the baryon density tl^h 2 and the matter 
density £l m h 2 . When the damping tail was first discovered by the CBI experiment ( Padin et al.|2001 j , 
it supplied compelling support for the standard theoretical modeling of the physics at recombination 
outlined here. Currently the best constraints on the damping tail are from the ACBAR experiment 
(Reichardt et al. 2008 see Fig. [7]). Consistency between the low order peaks and the damping tail 
can be used to make precision tests of recombination and any physics beyond the standard model at 
that epoch. For example, damping tail measurements can be used to constrain the evolution of the 
fine structure constant. 

Matter Power Spectrum: Shape & Amplitude — The acoustic peaks also determine the shape and 
amplitude of the matter power spectrum. Firstly, acoustic oscillations are shared by the baryons. In 
particular, the plasma motion kinematically produces enhancements of density near recombination 
(see Eqn.[IT3|) 

5b ~ —kr]*Vb{r]*) ~ — kr]*v-y(r]*) . (132) 
This enhancement then imprints into the matter power spectrum at an amplitude reduced by pi, / p m 



due to the small baryon fraction (Hu and Sugiyama 1996). Secondly, the gravitational potentials 



that the cold dark matter perturbations fall in are evolving through the plasma epoch due to the 
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Fig. 16. Transfer function. Acoustic and radiation physics is imprinted on the matter power spectrum as 
quantified by the transfer function T(k). The former is responsible for baryon oscillations in the spectrum and 
the latter a suppression of growth due to Jeans stability for scales smaller than the horizon at matter-radiation 
equality. 



processes described in §3.5| On scales above the horizon, relativistic stresses are never important and 
the gravitational potential remains constant. The Poisson equation then implies that 

A ~ (kr)) 2 <f>(0) , (133) 

and in particular, the density perturbation at horizon crossing where kr) ~ 1 is A = Ah ~ ^(0). For 
a fluctuation that crosses the horizon during radiation domination, the total density perturbation is 
Jeans stabilized until matter radiation equality 

^" 114 (w) Mpc - (134) 

Thereafter, relativistic stresses again become irrelevant and the potential remains unchanged until 
matter ceases to dominate the expansion 

$ « (^ eq )- 2 A H ~ (fc7? cq r 2 $(0) . (135) 

The transfer in shape from the initial conditions due to baryon oscillations and matter radiation 
equality is usually encapsulated into a transfer function T{k). Given an initial power spectrum of 



the form ( 94 ) , the evolution through to matter domination transforms the potential power spectrum 
to A; 3 P$/27r 2 oc k n ~ 1 T 2 (k) where T(k) oc k~ 2 beyond the wavenumber at matter-radiation equality. 
This scaling is slightly modified due to the logarithmic growth of dark matter fluctuations during the 
radiation epoch when the radiation density is Jeans stable. The matter power spectrum and potential 
power spectrum are related by the Poisson equation and so carry the same shape. Specifically 

where we have included a factor G(a) to account for the decay in the potential during the acceleration 
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epoch when relativistic stresses are again important (see e.g. Hu ( 2005[ )). This factor only depends 
on time and not scale as long as the scale in question is within the Jeans scale of the accelerating 
component. In this limit, G(a) is determined by the solution to 

d 2 G 



din a 2 



+ 4 + 



dlnH\ dG 



din a J din a 



+ 



3 + 



dlnH 
din a 



G = 0. 



(137) 



1 and G"(lna m d) = at an epoch a m d when the universe is 



with an initial conditions of G(lna m d) 
fully matter dominated. 

The transfer function T(k), with k in Mpc -1 , depends only on the baryon density Q^h 2 and 
the matter density Q m h 2 which are well determined by the CMB acoustic peaks. Features in the 
matter power spectrum, especially the baryon oscillations, then serve as standard rulers for distance 
measurements Eisenstein et al. (1998). For example, its measurement in a local redshift survey where 
the distance is calibrated in h Mpc -1 would give the Hubble constant h. Detection of the features 
requires a Gpc 3 of volume and so precise, purely local, measurements are not feasible. Nonetheless, 
the first detection of these oscillations by the SDSS LRG redshift survey out to z ~ 0.4 provide 
remarkably tight constraints on the acceleration of the expansion dEisenstein et al.||2005| ). 

The CMB also determines the initial normalization 5^ and so provides a means by which to test the 
effect of the acceleration on the growth function G(a) . The precision of this determination is largely set 
by reionization. The opacity provided by electrons after reionization suppress the observed amplitude 
of the peaks relative to the initial amplitude and hence 

4.6e~ (ai - T) x 10- 



(138) 



where r is the optical depth to recombination. Note that this is the normalization at k = 0.05 Mpc 1 
and even with uncertainties in the optical depth of 5t ~ 0.03 it exceeds the precision of the COBE 



normalization (c.f. Eqn. (95)). Finally, with the matter and baryon transfer effects determined, the 
acoustic spectrum also constrains the tilt. WMAP provided the first hints of a small deviation from 
scale invariance (Spergel et al. 2007) and the current constraints are n ~ 0.965 ± 0.015. 



Combining these factors into the conventional measure of the amplitude of matter fluctuations 
today 

dk k 3 P{k,a = 



4 



0"8 




0.563 



X (3.123/l) 



(139) 



where W a (x) 
8/i -1 Mpc. 



3x 3 (sinx — xcosx) is the Fourier transform of a top hat window of radius 



4 Polarization Anisotropy from Recombination 

Thomson scattering of quadrupolarly anisotropic but unpolarized radiation generates linear polar- 
ization. As we have seen in § |3.1[ £ > 2 anisotropy develops only in optically thin conditions. Given 
that polarization also requires scattering to be generated, the polarization anisotropy is generically 
much smaller than the temperature anisotropy. The main source of polarization from recombination 
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Fig. 17. _E-mode polarization power spectrum measurements. Adapted from BischofT et al. (20081 



is associated with the acoustic peaks in temperature. This source was first detected by the DASI 



experiment (Kovac et al. 20021 and recent years have seen increasingly precise measurements (see 



Fig. 17) 



Since acoustic polarization arises from linear scalar perturbations, they possess a symmetry that 
relates the direction of polarization to the wavevector or change in the polarization amplitude. As wee 



shall see in the next section, this symmetry is manifest in the the absence of i?-modes (Kamionkowski 



et al. 1997 Zaldarriaga and Seljak 1997). B- modes at recombination can be generated from the 



quadrupole moment of a gravitational wave. This yet-to-be-detected signal would be invaluable for 
early universe studies involving the inflationary origin of perturbations. 

We begin by reviewing the Stokes parameter description of polarization (£ |4.1[ ) and its relation 
to E and B harmonic representation (£ |4.2[ ). We continue with a discussion of polarized Thomson 
scattering in £4.3 In §4.4| and 4.5 we discuss the polarization signatures of acoustic oscillations and 



gravitational waves. 



4-1 Statistical Description 

The polarization field can be analyzed in a way very similar to the temperature field, save for one 
complication. In addition to its strength, polarization also has an orientation, depending on relative 
strength of two linear polarization states. 

The polarization field is defined locally in terms of Stokes parameters. In general, the polarization 
state of radiation in direction n described by the intensity matrix ^£ , j(n)£ , |(fi)^>, where E is the 
electric field vector in the transverse plane and the brackets denote time averaging. As a 2 x 2 
hermitian matrix, it can be decomposed into the Pauli basis 



P = c(E(n)Et(n)) 

= 9(n)cr + Q(n) cr 3 + U(h) cry + V(n) a 2 , 



(140) 
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where 

°"o = I n , I ■ <?2 





1 = i n > *3 = n n • (141) 



Orthogonality of the Pauli matrices says that the Stokes parameters are recovered from the polariza- 
tion matrix as Tr(<7iP) /2. The Stokes Q and f/ parameter define the linear polarization state whereas 
V defines the circular polarization state. We have chosen the proportionality constant so that all the 
Stokes parameters are in temperature fluctuation units. 

From this description, we see that Q represents polarization aligned with one of the principal axes 
of the transverse coordinate system whereas U represents polarization at 45° to these axes. 



4-2 EB Harmonic Description 

The disadvantage of the Stokes Q and U representation is that the distinction between the two 
depends on the coordinate system for the two transverse directions on the sky. Under a rotation of 
this basis by 9 

Q' ± ill' = e T2ie [Q ± ill] . (142) 

For a harmonic decomposition, it is more useful to choose this basis to be given by the wavevector 
itself. For small sections of the sky, the harmonic decomposition becomes a Fourier transform and 
we can define the E and B harmonics as 

E(\)±iB(\) = Jdh[Q'(h) ±iU'(h)]e- ilA (143) 

= e* 2i<tn J dh[Q(h) ± il/(n)]e- afi , 

where in the second line we have rotated Q and U back to a fixed coordinate system with the angle 
01 that the Fourier vector makes with the x axis. 

For linear scalar fluctuations, we have seen that the only direction for a given harmonic mode is 
set by the direction of the wavevector itself: velocity fields point in this direction or its opposite, 
quadrupole moments are symmetric about this axis, etc. This means that symmetry requires that 
such sources only generate Q' or E for each mode. This symmetry also holds once all the modes 
are superimposed back to the full polarization field: linear scalar perturbations generate only E- 
modes where the polarization direction is related to the direction in which the polarization amplitude 
changes. 

To generalize this decomposition to the full curved sky, we need to replace plane waves, the tensor 
eigenfunctions of the Laplace operator in a flat space, to the correct tensor eigenfunctions for the 
curved space. These are called the spin-2 spherical harmonics: 

V 2 ±2 Y £m [a 3 T iffi] = -[1(1 + 1) - 4] ±2 Y £m [<x 3 =p ia x \ . (144) 
They obey the usual orthogonality and completeness relations 

J dh s Y e * m (h) s Y em (h) = 5 u >5 mm , , 
]r s Y/ m (n) s Y £m (n') = 8{<j> - ^(cos 9 - cos 9') . (145) 
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Fig. 18. Thomson scattering geometry. A quadrupole anisotropy in the incoming radiation leads to linear 
polarization. For scattering at = 7r/2, only one component of the initially unpolarized radiation (E' 2 ) is 
scattered leaving one outgoing state (Ei) unpopulated. 



We can therefore decompose the linear polarization field just like the temperature field 

[Q(n) ± iU(h)] = - Y^i E em ± iBi m ]±2Yt m (n) . 

Urn 

Likewise the power spectra are given by 

and the cross correlation of E'-polarization with temperature by 
Other cross correlations vanish if parity is conserved. 



(146) 

(147) 
(148) 

(149) 



4-3 Thomson Scattering 

The differential cross section for Thomson scattering 

^ = A|E'.E|V, (150) 

is polarization dependent and hence scattering generates polarization. Here E' and E denote the 
incoming and outgoing directions of the electric field or polarization vector. Consider incoming 



radiation in the x direction scattered at right angles into the — y direction (see Fig. 18). Heuristically, 
incoming radiation shakes an electron in the direction of its electric field vector or polarization e' 
causing it to radiate with an outgoing polarization parallel to that direction. However since the 
outgoing polarization e must be orthogonal to the outgoing direction, incoming radiation that is 
polarized parallel to the outgoing direction cannot scatter leaving only one polarization state. 

The incoming radiation however comes from all angles. If it were completely isotropic in intensity, 
radiation coming along the z would provide the polarization state that is missing from that coming 
along x leaving the net outgoing radiation unpolarized. Only a quadrupole temperature anisotropy 
in the radiation generates a net linear polarization from Thomson scattering. As we have seen, a 
quadrupole can only be efficiently generated if the universe is optically thin to Thomson scattering a 
given perturbation. 
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Fig. 19. Polarized landscape. While the i?-spectrum and QE cross correlation are increasingly well measured, 
the .B-spectrum from inflationary gravitational waves (shown here near the maximal value allowed by the 
temperature spectrum) and gravitational lensing remains undetected. S hown here are pro j ected error bars 
associated with Planck sample variance and detector noise. Adapted from Hu and Dodelson (2002 1. 



4-4 Acoustic Polarization 

Acoustic oscillations in the dissipation regime provide the conditions necessary for polarization. Recall 
that radiative viscosity in the plasma is equivalent to quadrupole anisotropy in the photons. Since 



the quadrupole is of order (see Fig. 11 ) 



kv^ ( k \ w 7 

vr 7 ^ ~ — — , (151 

r \k D J koV* 

the polarization spectrum rises as I/Id to peak at the damping scale with an amplitude of about 10% 
of the temperature fluctuations before falling due to the elimination of the acoustic source itself due 
to damping. Since i> 7 is out of phase with the temperature, 

+ * oc cos(ks); u 7 oc sin(£;s) , (152) 

the polarization peaks are also out of phase with the temperature peaks. Furthermore, the phase 
relation also tells us that the polarization is correlated with the temperature perturbations. The cor- 
relation power Cj being the product of the two, exhibits oscillations at twice the acoustic frequency 

(0 + if!)(v~) oc cos(ks) sin(A;s) oc sin(2A;s) . (153) 



As in the case of the damping, the predicting the precise value requires numerical codes (Bond and 
Efstathiou|[l987 ) since f changes so rapidly near recombination. Nonetheless the detailed predictions 
bear these qualitative features. 

Like the damping scale, the acoustic polarization spectrum is uniquely predicted from the temper- 
ature spectrum once Q^h 2 , Sl m h 2 and the initial conditions are specified. Polarization thus represents 
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a sharp test on the assumptions of the recombination physics and power law curvature fluctuations in 
the initial conditions used in interpreting the temperature peaks. For example, features in the initial 
power spectrum appear more distinct in the polarization than the temperature due to projection 
effects (see e.g. Hu and Okamoto (2004)). 



4-5 Gravitational Waves 

During the break down of tight coupling that occurs at last scattering, any gravitational waves 
present will also imprint a local quadrupole anisotropy to the photons and hence a linear polarization 
to the CMB (Polnarev 1985). These contribute to the BB power and their detection would provide 
invaluable information on the origin of the fluctuations. Specifically, in simple inflationary models 
their amplitude gives the energy scale of inflation. The gravitational wave amplitude h oscillates and 
decays once inside the horizon, so the associated polarization source scales as h/f and so peaks at 
the I ~ 100 horizon scale and not the damping scale at recombination (see Fig. 19). This provides a 
useful scale separation of the various polarization effects. 

If the energy scale of inflation is near the 10 16 GeV scale then the signal is potentially detectable 
by the next generation of polarization experiments. 



5 Discussion 

The one and a half decades since the discovery of CMB anisotropy by COBE DMR has seen re- 
markable progress that ushered in the current epoch of precision cosmology. From the preliminary 
detections of degree scale power to the current measurements of 5 acoustic peaks, the damping tail and 
acoustic polarization, the milestones in the observational verification of our theoretical understanding 
of the universe at recombination have steadily been overtaken. Correspondingly, the measurements 
of the energy density contents of the universe at recombination and the distance to and hence ex- 
pansion rate since recombination have improved from order unity constraints to several percent level 
determinations . 

The next generation of experiments, including the Planck satellite and ground based polarization 
measurements, will push these determinations to the 1% level and beyond and perhaps detect the one 
outstanding prediction of the recombination epoch: the B-mode polarization of gravitational waves 
from inflation. Likewise the enhanced range of precision measurements of the power spectrum will 
bring measurements of the spectrum of scalar perturbations to the percent level and further test the 
physics of inflation. 

Beyond the recombination epoch reviewed here, these experiments will also test secondary tem- 
perature and polarization anisotropy from reionization, lensing, and galaxy clusters as well as the 
Gaussianity of the initial conditions. Combined, these measurements will test the standard cosmo- 
logical model with unprecedented precision. 

Acknowledgments: I thank the organizers, R. Rebelo and J. A. Rubino-Martin, as well as the students 
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